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I.  INTRODUCTION 

Future  applications  of  SHIVA-type  electromagnetically  driven  plasma 
Implosion,  Including  the  design  of  a  full-scale  nuclear  environmental  simu¬ 
lator,  will  require  the  storage  and  delivery  of  exceedingly  large  amounts  of 
electrical  energy.  As  presently  envisioned,  the  minimum  useful  radiation 
output  of  5-10  MJ  from  such  a  simulator  implies  delivered  electrical  energies 
of  15-30  MJ.  Using  design  philosophies  explored  in  previous  work  (Ref.  1), 
this  level  of  delivered  electrical  energy  implies  that  energies  of  30-60  MJ 
must  be  stored  in  a  form  from  which  they  can  be  delivered  in  microsecond  time 
scales.  Based  on  recent  estimates,  this  probably  results  in  a  system  costing 
$100  million  and  occupying  a  volume  of  54,000  m3  (70  m  dla  x  14  m  high). 

An  attractive  alternative  technology  for  supplying  large  amounts  of 
energy  may  be  conceptually  developed  using  inductive  (magnetic)  energy  con¬ 
ditioning  techniques  coupled  with  inertial  primary  energy  storage.  The 
inertial  primary  store  is  essentially  present  technology.  Therefore,  the 
purpose  of  this  report  is  to  explore  the  potential  applicability  of  inductive 
pulse-forming  techniques  for  driving  imploding  plasma  loads.  Special  atten¬ 
tion  is  paid  to  the  critical  problem  area  in  the  system,  which  is  the 
development  of  a  suitable  opening  switch  to  deliver  power  to  the  load. 

The  inductive  store  driving  a  cylindrical  implosion  system  is  particularly 
amenable  to  analytic  modeling,  and  the  first  section  of  this  report  discusses 
the  implication  of  such  models.  The  results,  while  encouraging  in  themselves, 
imply  rather  severe  operational  requirements  for  the  opening  switch.  These 
requirements  appear  well  beyond  currently  existing  technology.  The  second 
section  of  this  report  deals  with  the  investigation  of  more  detailed  models 
which  require  numerical  solutions.  The  more  detailed  models  address  more 
quantitatively  the  influence  of  switch  behavior  on  load  performance. 
Considerably  more  latitude  is  allowed  than  is  suggested  by  simple  analytical 
models.  The  last  section  discusses  existing  data  on  the  performance  of 
opening  switches  in  light  of  the  performance  required  by  the  models. 


1.  Reinovsky,  R.  E.,  Conceptual  Design  of  a  6  MJ  Pulse  Power  System  for 
Driving  SHIVA,  AFWL-DYP-TN-78-126,  Air  Force  Weapons  Laboratory, 
Kirtland  Air  Force  Base,  New  Mex,  December  1978. 
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II.  ANALYSIS 

This  section  explores  the  concept  of  inductively  driven  imploding  loads 
using  the  simplest  analytic  circuit  models.  This  analysis  differs  slightly 
from  the  familiar  analysis  of  the  inductive- store/ inductive-load  problem 
in  that  it  allows  for  a  load  whose  inductance  increases  subsequent  to  the 
initiation  of  current  interruption  in  the  primary  circuit.  The  practical 
circuit  under  consideration  is  shown  in  Figure  1.  It  consists  of  a  d.c. 
charged  capacitor  bank  (C)  discharged  through  an  inductor  (L$)  and  a  switch 
which  is  initially  closed  but  will  open  to  transfer  energy.  The  inductor 
includes  the  internal  inductance  of  the  bank  (Li )  and  may  be  augmented  by 
additional  storage  inductance  (Lgxt).  Thus,  the  energy  stored  in  the  para¬ 
sitic  bank  inductance  is  made  available  to  the  circuit.  An  initially  open 
output  switch  S2  isolates  the  load  from  the  charging  system  prior  to  switching 
time  t$.  The  load  inductance  consists  of  a  nonzero  initial  inductance 
Lq  and  a  time  varying  inductance  AL^  corresponding  to  the  implosion  of  the 
foil /plasma.  The  aL  in  nH  is  given  by 

aL(tpj  • 2  h  ,n  (^) 

where  is  the  final  radius  of  the  pinch,  RQ  is  the  initial  radius  of  the 
foil  cylinder,  and  h  is  the  height  of  the  foil  cylinder  in  cm.  The  ratio 
RQ/Rf,  referred  to  as  the  convergence  ratio,  is  controlled  by  magnetohydro¬ 
dynamics  (MHD)  considerations  and  typically  lies  between  10  and  20.  For 
implosion  times,  greater  than  1  ys  a  convergence  of  7-12  is  common;  and,  for 
shorter  time  scales,  improved  stability  is  expected  to  lead  to  slightly 
higher  convergences. 

In  operation,  energy  is  transferred  from  C  to  L$  and  Sj  is  opened  at  peak 
current  (IQ).  At  peak  current  the  voltage  on  C  is  near  zero;  and,  if  the 
energy  transfer  to  the  load  is  fast  compared  to  the  LSC  time  constant,  then 
the  amount  by  which  the  charge  on  C  can  charge  (i.e.,  AQ  =  IQAt)  is  small  and 
the  voltage  on  C  can  be  assumed  to  be  near  zero  and  to  remain  near  zero  during 
all  time  of  interest.  Thus,  the  practical  circuit  in  Figure  1  is  reduced  to 
the  equivalent  circuit  in  Figure  2. 
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Figure  1.  General  inductive  energy  storage  circuit. 


Figure  2.  Equivalent  circuit. 


Approach  the  analysis  of  the  circuit  in  Figure  2  in  a  two-step  process. 
In  each  step,  conservation  of  magnetic  flux  is  assumed.  The  first  step  is 
the  opening  of  the  switch  Si  which  is  considered  a  resistor  whose  resistance 
goes  from  zero  to  infinity  rapidly.  The  second  step  is  the  implosion,  with 
Si  effectively  out  of  the  circuit  and  S2  closed,  in  which  goes  from  Lq 
to  (Lq  +  aL).  Prior  to  the  opening  of  Si  the  energy  stored  in  L$  is 


Eo  = 


TV 


(1) 


After  Si  opens  and  S2  closes  at  t$,  the  current  flowing  assumes  a  new  value 
Ii  given  by 


Ii  = 


Ls  + 


(2) 
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Ii  is  flowing  in  both  L$  and  where  L(ts)  =  LQ.  Thus  the  energy  stored  in 
L 

L. 


inductors  L  and  is 
s  o 


E,  = 


1  '  Ls  +  Lo  0 


(3) 


and  the  energy  dissipated  in  the  switch  resistance  as  it  climbs  to  infinity  is 

L 


E 


(4) 


switch  l.  +  L.  o 
s  o 

Note  that  to  minimize  the  energy  lost  in  the  interruption  and  transfer  LQ 
should  be  as  low  as  possible. 

After  the  implosion,  the  value  of  L ^  has  risen  to  Ln  +  aL,  the  current 


has  dropped  to  I2, 


h  = 


2  “  L  +  L  +  aL  lo 

S  0 

and  the  energy  remaining  in  the  magnetic  field  of  L$,  Lq  and  aL  is  only 


(5) 


E* '  l."¥  rva.  Eo  (6> 

The  difference  in  energy  E:  -  E2  is  the  energy  dissipated  in  doing  work  on 
the  load  inductance  to  increase  its  value  from  LQ  to  (LQ  +  aL).  Thus, 

Ex  -  E2  is  just  the  kinetic  energy  coupled  to  the  plasma  shell  during  the 
implosion. 


L*  =  '  Ls*Lot4L) 


(7) 


A  kinetic  efficiency  can  be  defined  as  the  ratio  of  kinetic  energy  to 
inductively  stored  energy.  Thus 


Eke  .  Ls 


'ke  E. 


L+L  L+L+aL 
s  o  s  o 


(8a) 


_ _ 

'ke  ‘  (L$  +  Lq)  (L$  +  L0  +  AL) 


(8b) 
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Recalling  that  LQ  «  L$  is  desired,  LQ  compared  to  L$  can  be  neglected,  and 


aL 

nke  *  Ls  +  AL 


(9) 


Thus,  if  aL  can  be  made  much  larger  than  L$  (Ls  <<  aL),  then  nke  -*■  1,  and  we 
can  in  principle  get  virtually  all  the  magnetically  stored  energy  into  kinetic 
energy.  Thus,  in  general,  high  efficiency  requires 

Lq  «  L$  «  AL  (10) 


Lq  is  in  general  fixed  by  the  mechanical  and  electrical  consideration  of 
power  flow  to  the  load,  while  AL  is  fixed  by  the  convergence  ratio  and  the 
length.  Usually  values  of  LQ  below  3-5  nH  are  unrealistic,  while  conversion 
ratios  of  20  and  a  length  of  a  few  (1-4)  cm  results  in  AL  of  6-24  nH.  Thus, 
in  practice,  the  criteria  cannot  be  easily  satisfied  to  produce  very  high 
efficiency. 

The  value  of  l$  remains  as  one  parameter  which  can  be  adjusted.  If  L$ 
is  chosen  to  be  very  small,  a  large  amount  of  energy  is  lost  in  the  switching 
operation;  if  it  is  chosen  to  be  large,  the  energy  transfer  to  the  load 
suffers.  To  find  the  optimum  choice  of  L$  we  take  dnke/dLs  and  set  the  result 
equal  to  zero.  This  results  in  a  criterion  for  L$,  namely: 

Ls  ■  vV  *  Lo  4L  Cl) 


For  the  case  of  the  static  load,  aL  =  0,  Equation  11  gives  Lg  =  LQ.  Thus 
the  familiar  static  result  is  recovered;  and,  as  expected  from  Equations  2-4, 


Ei  =  50%  Eq 


(12a) 


Eswi tch 


=  50%  Eq 


(12b) 


Figure  3  is  a  plot  of  nke  as  a  function  of  Lg  for  an  implosion  where 
aL  =  12  nH  and  for  values  of  L  of  1,  3,  and  5  nH.  The  figure  shows  that  the 
efficiency  of  coupling  inductive  to  kinetic  energy  goes  through  a  maximum  at 
the  predicted  optimum  L,.  and  that  the  optimum  is  quite  broad,  generally  losing 
less  than  10  percent  of  peak  efficiency  for  a  factor  of  2  deviation  from 
optimum  storage  inductance.  Thus,  while  designing  to  the  optimum  inductance 
is  useful,  precise  matching  is  usually  unnecessary. 
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Equations  2-4  show  that  achieving  values  of  LQ  smaller  than  L$  Is 
desirable,  but  that  significantly  reducing  LQ  Is  one  of  the  harder  tasks  In 
system  designs.  The  gain  in  efficiency  to  be  achieved  from  lowering  LQ  must 
be  evaluated  quantitatively.  Figure  4  is  a  plot  of  kinetic  efficiency  from 
Equation  8  (kinetic  energy  per  unit  Inductively  stored  energy)  versus  AL  for 
a  range  of  LQ  from  1-10  nH.  As  expected,  nke  increases  with  increasing  aL 
and  decreasing  LQ.  For  aL  of  12  nH  (2  cm  high,  convergence  of  20)  the  values 
of  L  *  1  give  a  very  respectable  efficiency  of  about  57  percent,  while  for 
Lq  =  10  nH,  the  efficiency  drops  to  less  than  20  percent.  This  Indicates 
that  a  decrease  in  LQ  is  an  effective  method  of  Increasing  efficiency. 

The  other  method  of  Increasing  nkg  is  by  increasing  aL.  Since  the  con¬ 
version  ratio  is  essentially  fixed,  aL  can  be  increased  by  increasing  the 
height,  generally  at  little  or  no  cost  in  increasing  LQ.  For  example, 
doubling  the  height  of  a  2-cm  load  to  4  cm  doubles  aL  from  12-24  nH,  which 
for  Lq  *  5  nH  Increases  nke  from  30  to  42  percent.  For  those  applications 
where  kinetic  energy  per  unit  length  is  a  criterion  in  addition  to  total 
kinetic  energy,  the  process  of  lengthening  the  foil  actually  results  in  a 
decrease  in  kinetic  energy  per  unit  length,  and  is  therefore  not  a  useful 
approach. 

For  electrically  exploded  conductors,  the  operation  of  the  switch  is 
determined  by  the  energy  (and  to  some  extent  power  history)  dissipated  in  the 
switch.  Therefore  it  is  useful  to  characterize  the  circuit  performance  in 
terms  of  the  energy  dissipated  in  the  switch.  Equation  4  gives  Eswitch- 
Figure  5  is  a  plot  of  the  dissipation  fraction  5  defined  as  the  ratio  of 
energy  dissipated  in  the  switch  to  the  total  inductively  stored  energy 


Equation  13  indicates  that  6  is  a  function  of  aL  through  L$  when  L$  is 
optimized  by  Equation  11,  and  is  therefore  a  relatively  weak  function  of  Al. 
For  small  AL  the  curve  approaches  50  percent,  as  expected  for  a  static  load. 
The  significance  of  5  is  that  it  represents  the  minimum  amount  of  energy 
that  will  be  dissipated  when  the  switch  opens,  regardless  of  the  character¬ 
istics  or  relative  time  of  the  operations  of  switch  Sj  and  switch  S2. 
Conversely,  it  is  the  minimum  energy  available  to  actuate  a  dissipation 
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Figure  4.  Kinetic  efficiency  vs  change  in  load  inductance. 
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driven  switch.  There  is  a  temptation  to  develop  a  switch  which  requires  very 
little  energy.  Figure  5  indicates  that,  for  example,  for  Lfl  *  3  nH  and 
aL  -  12  nH,  almost  40  percent  of  the  energy  must  go  into  the  switch,  regard¬ 
less  of  how  clever  the  design. 

Since  dissipatively  operated  current  interrupting  switches  are  likely 
candidates,  and  since  a  dissipative  switch  may  require  more  than  the  minimum 
energy  <5,  the  fraction  f  of  the  inductively  stored  energy  coupled  to  kinetic 
energy  after  switching  is  also  a  relevant  parameter.  Using  Equations  2,  3, 
and  7 


F  = 


o _ 

+  AL 


(14) 


Figure  6  is  a  plot  of  the  coupling  fraction  f.  The  plot  shows  that,  for 
realistic  values  of  3  nH  <  LQ  <  5  nH  for  aL  of  12  nH,  approximately  half  the 
energy  remaining  after  switching  is  coupled  to  kinetic  energy. 

When  Si  opens  and  S2  closes,  the  current  flowing  in  the  charging  circuit 
falls  from  I  to  the  value  given  by  Equation  2. 


The  current  in  the  load  rises  from  0  to  1^.  The  rate  at  which  the  load 
current  rises  may  be  important  from  the  standpoint  of  MHD  stability  of  the 
imploding  plasma.  If  the  switch  resistance  rises  from  R$  =  0  to  R$  s  R 
instantly  (a  step  function),  and  if  a  form  1^  =  I,  (l-e"ut)  is  assumed  for 
IL,  then  we  observe  from  Figure  2  that  the  voltage  across  the  resistor 
(a  short  time  At  following  the  step  increase  in  resistance)  is 

VR  =  roR  =  LL  IT  ( 


Thus 


A  loR 

dt 


=  ali  e 


■at 


For  e 


■at 


!o  R 
“  =  Ii  Li 


(16) 


(17) 
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or  by  Equation  2, 


The  time  constant  reflected  In  Equations  17  and  18  may  also  be  found  by 
simultaneously  rating  the  two  loop  equations  for  the  current  in  Figure  2. 

That  solution  verifies  both  the  assumed  form  of  I,  and  the  value  of  a.  For 
very  large  values  of  R,  the  rise  time  in  the  load  is  very  short.  For  example, 
if  R  *  1  fl  and  LQ  *  5  nH,  Ls  =  10  nH,  a  =  0.3  x  109,  the  current  rise  time 
would  be  3  ns. 

From  this  simple  analysis,  design  criteria  emerge: 

a.  Minimize  L  . 

o 

b.  Choose  Lc  «  \l L  *  +  L  aL  . 

s  v  o  o 

c.  Maximize  aL,  but  not  by  increasing  h  if  kinetic  per  unit  length  is 
a  concern. 

d.  A  significant  amount  (at  least  25  percent  in  most  cases)  of  energy 
will  be  dissipated  in  the  switch,  placing  a  lower  limit  on  the  switches 
dissipation  rating. 

This  analysis  makes  a  number  of  assumptions: 

a.  Infinite  final  resistance  of  the  opening  switch  is  assumed.  In  a 
crude  attempt  to  quantify  "infinite,"  the  final  resistance  of  the  switch 
must  be  large  compared  to  the  highest  impedance  presented  by  the  load  to 
prevent  significant  current  sharing  through  the  switch.  The  L  component 
given  by 


i  _  °  h 
L(t)  '  ST  h 


(19) 


(where  v  and  r  are  velocity  and  position,  respectively,  of  the  foil)  is  the 
largest  contributor  to  the  load  impedance.  Thus,  achieving  a  final  switch 
resistance  much  greater  than  Lisa  possible  criterion. 


b.  It  was  assumed  that  the  resistance  rises  quickly.  Quantitatively 
quickly  implies  that 


Atswitch  <<  At1mplosion 


(20) 
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c.  It  is  assumed  that  the  output  switch  S2  is  closed  when  the  resistance 
of  Si  begins  to  rise.  Since  total  time  of  current  flow  in  the  load  may  be  a 
determiner  of  foil /plasma  stability,  It  may  be  desirable  to  close  the  switch 
later,  thereby  achieving  higher  initial  voltages  at  closure  time  and  hence 
presumably  a  higher  initial  i  In  the  load. 
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III.  CIRCUIT  CALCULATIONS 

The  assumptions  outlined  In  Section  II  place  severe  constraints  on  system 
components,  especially  on  the  opening  switch.  The  validity  of  the  require¬ 
ments  placed  by  the  constant  flux  model  can  be  evaluated  using  a  more  detailed 
circuit  model.  In  this  section,  the  numerical  solution  to  a  circuit  similar 
to  that  In  Figure  1  will  be  discussed.  For  numerical  solutions,  values  were 
chosen  for  circuit  elements  which  correspond  to  the  parameters  of  the  SHIVA-I 
capacitor  bank  systems  as  shown  in  Figure  7. 


Figure  7.  Circuit  used  for  numerical  circuit  solutions. 

The  capacitor  is  267  uF,  charged  to  120  kV  storing  1.9  MJ.  The  load  is 
modeled  as  a  time  changing  inductance  with  values 

LLD  =  ^  h  In 


where  R  is  the  radius  of  the  return  conductor  and  r ^  is  the  initial  radius 
of  the  imploding  foil;  and  the  time  varying  resistance  has  a  value  of 


RLD  =  ip-  h  (22 

2*  r(t) 

where  Is  the  velocity  of  the  imploding  cylinder.  The  coupled  differen¬ 
tial  equations  of  motion 
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complete  the  set  that  Is  to  be  solved  numerically.  The  Initial  foil  radius 
was  chosen  at  5  cm  and  the  height  at  2  cm  for  stability  reasons.  The  return 
conductor  radius  was  chosen  at  17.5  cm  to  give  an  Initial  value  for  the  LLD 
of  5  nH,  which  correspond  to  LQ  =  5  nH  in  the  analytic  model.  The  assumption 
of  20:1  conversion  leads  to  a  minimum  radius  of  2.5  mm,  a  final  value  LLD  of 
17  nH,  and  a  a L  of  12  nH.  Thus  from  LQ  and  aL  a  value  of  the  storage  induct¬ 
ance  is  chosen  from  Equation  11  to  be  9.2  nH.  The  series  output  switch  is 
modeled  as  R2,  which  is  given  as  a  time  varying  resistor  whose  value  is  1  Mn 
prior  to  switching  time  Tg  and  changes  to  0.1  mn  in  5  ns  subsequent  to  T  . 

The  current  interrupting  switch  is  Ri,  whose  value  is  varied  as  a  problem 
variable.  The  resistance  of  Ri  is  arbitrary  and  is  not  computed  from  a 
material  model  for  this  analysis.  The  circuit  was  subjected  to  numerical 
analysis  using  a  circuit  solving  code  called  Sceptre  for  a  variety  of  Ri 
profiles  and  time  scales  and  for  a  variety  of  switch  times  T  . 

TERMINAL  RESISTANCE 


For  the  circuit  in  Figure  7  with  LQ  =  5  nH,  and  AL  =  12  nH,  the 
efficiency  from  Figure  4  is  30  percent.  Assuming  all  capacitively  stored 
energy  (1.9  MJ)  is  transferred  to  the  inductance  Lg,  the  coupled  kinetic 
energy  should  be  570  kJ.  With  a  desired  final  velocity  of  33  cm/ys,  the 
total  foil  mass  should  be  approximately  1  x  10“ 5  kg,  or  161  ug/cm2.  At  33 
cm/ys  final  velocity  and  2.5  mm  final  radius,  the  final  L  is  0.5  n.  The 
implication  of  the  constant  flux  model  was  that  the  final  valued  Rj  should 
be  much  greater  than  0.5  n  to  assure  that  most  of  the  current  is  flowing  in 
the  load.  To  model  the  situation,  the  resistance  profile  shown  in  Figure  8 
was  Imposed  on  Rj . 
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Where  RQ  was  taken  as  0.0  and  Ti  was  taken  as  the  time  of  peak  current, 

Ti  =  i  VT^T  3  2.46  ps  (24) 

To  assure  that  At  «  t^mp,  At  was  taken  as  100  ns  and  R^  was  varied  from 
30  mn  to  5  fl.  The  time  at  which  the  output  switch  closed  Tg,  was  a  constant 
at  2.465  ps.  Figure  9  shows  a  plot  of  the  kinetic  energy  coupled  to  the 
imploding  foil  and  the  final  velocity  of  the  foil  when  it  had  collapsed  to 
a  radius  of  2.5  mm  as  a  function  of  R^,  the  final  switch  resistance.  As 
anticipated,  the  kinetic  energy  coupled  at  lower  values  of  R^  is  lower  than 
that  observed  at  higher  values.  Somewhat  surprising  Is  the  fact  that  when 
Rf  «  n/2  (peak  L),  nearly  90  percent  of  the  kinetic  energy  predicted  by  the 
flux  model  was  observed  coupled  in  the  numerical  solution;  but  when  Rf  drops 
more  than  an  order  of  magnitude  to  30  mn  the  kinetic  efficiency  decreases 
only  moderately  to  62  percent  of  the  efficiency  predicted  by  the  flux  model. 
This  relatively  moderate  impact  of  changing  Rf  can  be  motivated  by  referring 
to  Figure  10  where  RLD  and  Ri(t)  are  plotted  as  functions  of  time  for  the 
case  where  Rf  3  30  mn.  The  plot  shows  that  the  dissipative  part  of  the  load 
impedance  RLD  rises  rapidly  at  the  very  end  of  the  implosion,  and  even  for  a 
very  modest  value  for  Rf,  RLD  is  less  than  Rf  for  about  90  percent  the  implo¬ 
sion  time.  Although  most  of  the  kinetic  energy  is  coupled  late  in  the 
implosion,  note  that  once  transferred  to  the  inductance  LLD,  energy  will  tend 
to  dissipate  more  rapidly  in  the  larger  resistance  (RLD)  because  of  the 
shorter  L/R  times.  In  other  words,  the  energy  lost  in  Rx  is  I2  Ri  At  (where 
At  Is  very  short),  and  thus  can  account  for  only  a  small  fraction  of  the 
inductively  stored  energy.  Thus,  since  RLD  is  much  less  than  Rf  for  most  of 
the  implosion,  efficient  coupling  to  LLD  is  achieved,  and  most  of  that 
energy  is  still  made  available  to  the  implosion.  Conceptually  carrying  this 
logic  to  the  extreme,  the  switch  operation  could  be  carried  out,  energy 
stored  in  LQ,  and  Ri  could  even  be  crowbarred  while  the  foil  implodes.  A 
lower  limit  on  kinetic  efficiency  can  be  found  by  considering  Equation  3 
for  the  total  energy  after  switching  Ei  =  (L$/Ls  +  LQ)  EQ.  The  fraction 

stored  in  L„  is 
o 


s  Lo 


(25) 
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Figure  10.  Load  L  vs  time. 


And  the  energy  coupled  to  kinetic  energy  is 

Lo  \  -  aL  _  LsLo 

'Lo  +  aL/ \  "  <Lo  +  aL>  (LS  +  Lo)Z  0 

LsLo  AL 

nCB  =  (Lq  +  AL)  (L$  +  Lq)2 

For  LQ  =  5  nH,  aL  =  12  nH,  and  L$  =  9.2  nH, 
riQg  =  16  percent 


(26) 

(27) 


(28) 


For  =  10  mn.  Figure  9  shows  nke  is  even  less  than  the  crowbarred  efficiency. 
This  is  explained  by  realizing  that  in  the  case  of  very  low  Rx ,  current 
transfer  to  the  load  is  never  complete.  The  implication  is  that  if  high  final 
fuse  resistance  is  not  easily  achievable,  a  crowbarred  case  when  R^.  goes 
temporarily  high  and  then  collapses  may  be  preferable  to  a  lower  fuse 
resistance. 
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Figure  9  also  exhibits  a  falloff  of  nke  above  approximately  0.5  n.  This 
result,  perhaps  more  surprising  than  the  relatively  moderate  falloff  at  low 
Rf,  is  explained  by  referring  to  Figure  11  where  the  energy  dissipation  in 
the  fuse  at  T  =  2.56  us,  the  time  when  the  resistance  stops  rising,  and  the 
energy  dissipated  in  the  fuse  at  pinch  time  are  plotted  as  a  function  of  Rf. 
The  simple  analysis  indicates  that  670  kJ  is  the  minimum  energy  that  must  be 
dissipated  in  the  fuse.  For  values  of  Rf  between  100  irft  and  1  n,  the  circuit 
puts  most  of  the  dissipated  energy  into  the  fuse  during  the  opening  time,  and 
little  energy  into  it  during  the  implosion.  This  behavior  is  similar  to  that 
of  the  simple  analysis.  This  range  also  corresponds  to  the  region  of  peak 
kinetic  energy  shown  in  Figure  9. 


At  larger  values  of  R^,  Figure  11  shows  that  excessive  energy  is  dissipated 
in  the  fuse  during  opening  time,  thus  leaving  less  energy  in  the  magnetic 
circuit  to  drive  the  implosion;  this  explains  the  lower  overall  efficiency. 

The  excessive  energy  deposited  at  high  values  of  Rf  was  investigated. 

Decreasing  the  closing  time  of  the  series  output  switch  R2  resulted  in  improved 
efficiency  and  decreased  dissipation  at  large  values  of  Rf.  In  fact,  the 
kinetic  energy  approaches  the  670  kJ  flux  model  value.  Presumably  when  R  is 
very  large  (i.e.,  R^  is  large  and  At  is  fixed),  excessive  energy  is  dissipated 
in  Ri  prior  to  the  closing  of  R2.  Actually  the  5  ns  value  of  R2  used  to 
generate  the  data  in  Figure  9  is  more  representative  of  practical  multichannel 
switches  than  the  <  1  ns  value  required  to  achieve  flux  model  efficiency.  The 
implication  is  that  low  jitter  output  switching  is  required  if  very  large 
values  of  Rf  are  used.  Moderate  values  of  R^  do  not  place  severe  limitations 
on  R2  performance.  The  conclusion  is  that,  for  implosion  parameters  discussed, 
values  of  Rf  that  are  greater  than  the  initial  L  of  the  load  (a  few  milliohms), 
but  not  much  greater  than  the  final  L  (one-half  ohm)  performance,  seem  to  be 
predicted  by  the  simple  model  within  about  20  percent.  Thus  the  criteria 
result 


Lo  » 


Rf. 


■pinch 


(29) 


OPENING  TIME 

The  simple  flux  model  presumes  that  the  implosion  is  carried  out  in  two 
steps.  First  a  current  interruption  occurs;  then  an  implosion  phase  occurs. 
The  energy  transfer  is  calculated  on  the  assumption  that  does  not  change 
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during  the  interruption  phase;  that  is,  the  calculation  assumes  switching  into 
a  static  load.  Figure  12  is  a  plot  of  LLD  vs  time  for  a  case  where  R^.  =  100 

mfi,  aL  =  12  nH,  LQ  =  5nH,  L$  =  9.2  nH,  and  AL  =  100  ns.  During  the  300  ns 

following  output  with  closure  at  2.465  ys  (2.465  to  2.765  us),  the  inductance 
of  the  load  LLD  rises  about  25  percent.  Thus,  for  time  scales  of  about  three- 
fourths  of  the  implosion  time,  the  opening  switch  is  seeing  a  nearly  constant 
L^  approximately  equal  to  L^t  _  g^.  Figure  13  is  a  plot  of  kinetic  energy 
coupled  for  a  20  to  1  compression  vs  the  ramp  time  of  the  switch  impedance. 

The  final  value  R^  is  100  mn.  The  plot  shows  virtually  no  change  of  efficiency 

for  time  scales  up  to  300  ns,  which  is  very  close  to  75  percent  of  the  implo¬ 

sion  time.  For  time  scales  up  to  500  ns,  the  loss  of  efficiency  is  less  than 
5  percent  and  the  implosion  time  lengthens  somewhat  (from  450-550  ns).  For 
ramp  times  (At)  longer  than  500  ns,  the  switch  can  be  thought  of  simply  as  a 
500-ns  ramp  of  lower  R  and  thus  a  lower  R^.  The  effect  of  lowering  Rf  is 
shown  in  Figure  9. 

IMPLOSION  MASS/FINAL  VELOCITY 

One  of  the  advantages  of  inductively  driven  implosion  systems  is  that,  at 
least  in  the  simplest  model,  the  kinetic  energy  coupled  is  dictated  only  by 
the  inductance  ratios  and  is  independent  of  the  implosion  mass.  This  allows 
relatively  wide  variations  in  final  velocity  to  be  achieved  independently  of 
kinetic  energy,  and  hence  allows  assessment  of  the  effect  of  final  velocity 
on  the  thermal izati on  process.  Figure  14  is  a  plot  of  the  kinetic  efficiency, 
final  velocity,  and  implosion  time  as  a  function  of  implosion  mass.  The  plot 
shows  that,  for  a  full  order  of  magnitude  change  of  implosion  mass  (5  x  10“5 
to  5  x  10"6  kg),  the  change  in  kinetic  energy  is  small  while  the  change  in 
velocity  is  given  by  the  anticipated  vTO  from  11.9  cm/us.  As  expected,  the 
implosion  time  varies  over  a  similarly  wide  range  associated  with  the  changing 
final  velocity.  For  very  large  masses  and  for  very  small  masses,  the  kinetic 
efficiency  suffers  somewhat.  The  numerical  solutions  show  that  for  the  large 
masses  the  long  implosion  time  leads  to  reverse  charging  of  the  bank  capaci¬ 
tance  (because  a  relatively  large  current  is  flowing  in  the  positive 
direction  for  a  long  time  after  current  peak).  The  energy  stored  in  the 
recharging  capacitor  is  approximately  three  times  the  observed  loss  in  kinetic 
energy.  For  small  values  of  mass  (hence  short  implosion  time),  loss  results 
from  effective  lowering  of  the  fuse  resistance  caused  by  early  pinch. 
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Figure  12.  Load  L  vs  time. 
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Figure  13.  Kinetic  energy  vs  rise  time. 
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OUTPUT  SWITCH  CLOSURE  TIME 

The  time  at  which  the  output  switch  closes  is  expected  to  impact  the  total 
kinetic  efficiency.  For  example,  closing  the  output  switch  late  in  the  inter¬ 
ruption  may  be  expected  to  result  in  excessive  energy  deposition  in  the  fuse 
and  hence  lower  kinetic  efficiency.  On  the  other  hand,  closure  of  the  output 
switch  too  early  may  be  expected  to  result  in  lower  voltages  across  the  load 
(and  hence  lower  initial  I)  and  perhaps  a  longer  implosion  time  for  a  given 
load. 

Figure  15  shows  a  plot  of  nke  and  implosion  time  as  function  switch  time 
T$  for  a  case  where  R^  equals  500  mfi.  The  implosion  mass  was  1  x  10“5  kg, 
and  the  switch  opening  time  At  was  100  ns.  For  reference,  the  fuse  resistance 
profile  is  also  sketched.  The  figure  shows  that  both  kinetic  energy  and 
implosion  time  are  sensitive  to  switch  closure  time.  As  expected,  kinetic 
energy  drops  and  implosion  time  increases  with  later  closing  time.  Implosion 
time  shows  a  tendency  to  flatten  out  for  closure  time  near  the  start  of  the 
interruption  (2.46  ys).  Fortunately,  from  a  practical  point  of  view,  the 
earliest  possible  closure  time  (after  start  of  interruption  time)  appears 
most  promising  from  both  the  point  of  view  of  efficiency  and  of  implosion  time. 
This  suggests  that  a  self-breaking  switch  set  to  close  at  low  (perhaps  a  few 
x  104  V)  may  be  adequate.  The  difficulty  is  that  of  achieving  low  inductance 
multichannel  performance  with  low- voltage  self -breaking  switching. 

INITIAL  FUSE  RESISTANCE 

Two  cases  were  compared  with  R^  *  100  mo;  one  with  Ri  =  0  and  one  with  R. 

■  1  mo.  Significantly  less  (300  kJ)  kinetic  energy  was  observed  with  R^ 

~  1  mfi  then  with  R^  =  0.  On  the  other  hand,  peak  storage  current  was  also 
reduced,  as  expected  with  the  larger  values  of  R^;  but  the  kinetic  efficiency 
was  virtually  identical  in  the  two  cases.  The  conclusion  is  that  kinetic 
efficiency  is  virtually  independent  of  changes  in  the  ratio  Rf/R.j . 

RESISTANCE  SHAPE 

The  linear  ramp  is  conceptually  the  easiest  to  examine,  but  is  not 
necessarily  the  most  accurate  representation  of  resistivity  data.  A  number 
of  profiles  were  compared,  each  of  which  rose  from  =  0.1  mo  to  500  mfi 
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in  At  *  200  ns.  The  ramp  profile  was  compared  with  a  smoother  S-shaped  curve, 
a  square  low  profile,  and  an  exponential.  All  four  profiles  produced  from 
488-523  kJ  coupled  kinetic  energy  or  a  7  percent  variation.  The  linear  ramp 
produced  less  kinetic  energy  than  either  the  S  shape  or  the  exponential,  and 
slightly  more  than  the  square  low  profile. 
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IV.  SWITCH  TECHNOLOGY 

Finally,  it  is  appropriate  to  consider  the  prospects  for  the  success  of  a 
high-energy  inductive  store/opening  switch  system  as  a  driver  for  a  practical 
imploding  plasma  load.  Significant  data  have  been  published  on  the  behavior 
of  exploded  foil  fuses  used  as  opening  switches;  but  in  general  the  energy 
level  (25  kJ)  and  the  time  scale  (10  to  a  few  hundred  us)  are  not  representa¬ 
tive  of  the  behavior  of  the  fusing  element  in  systems  of  interest  (2  MJ,  1-2 
us).  The  work  most  nearly  approaching  these  parameters  is  that  already 
performed  by  the  Air  Force  Weapons  Laboratory  (AFWL)  at  the  100-kJ,  3-4-s 
level  (Ref.  2).  This  section  examines  the  results  of  this  work  In  light  of 
the  foregoing  analysis  and  circuit  calculations. 

Figures  16  and  17,  extracted  from  Reference  2,  show  current  and  voltage 
profiles  for  a  set  of  copper  foil  fuses  quenched  in  glass  beads  for  a  variety 
of  physical  lengths  and  widths  which  maintain  a  constant  total  fuse  mass  of 
25  g  (for  1  mil  thickness).  Based  on  preceding  analysis,  the  most  promising 
choice  for  a  fuse  might  be  one  which  produces  the  highest  storage  current 
while  still  opening  in  times  less  than  (but  not  necessarily  much  less  than) 
the  implosion  time.  From  the  data,  it  is  clear  that  (for  this  family) 
opening  time  must  be  compromised  against  peak  current.  It  is  also  clear  that 
opening  time  does  not  have  a  simple  interpretation.  It  is  convenient  to 
accept  the  FWHM  of  the  voltage  pulse  to  measure  quality  of  the  switch. 
Consideration  of  the  simple  L-R  circuit  shows  that  the  FWHM  of  the  voltage 
pulse  is  composed  of  two  parts.  The  rise  time  is  related  to  the  rise  time 
of  the  fuse  resistance.  In  the  simplest  case  of  a  step  function  resistance, 
the  voltage  rise  is  also  a  step;  however,  the  fall  time  is  related  to  the  L/R 
decay  of  the  circuit  L  and  the  final  switch  resistance  R.  Thus,  a  sharp 
rising  pulse  with  a  fast  decay  indicates  both  rapid  dR/dt  and  high  final 
resistance,  hence  a  better  quality  switch. 


2.  McClenahan,  et  al.,  200  Kilojoule  Copper  Foil  Fuses,  AFWL -TR- 78- 130, 
Air  Force  Weapons  Laboratory,  Kirtland  Air  Force  Base,  New  Mex, 

(to  be  published). 
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Figure  15  shows  that  the  shortest  interrupt  time  is  associated  with  the 
peak  voltage  (maximum  I)  but  not  with  the  maximum  storage  current;  thus,  com¬ 
promise  is  in  order.  For  this  analysis,  the  maximum  voltage  case  Is  discussed. 
The  FWHM  of  this  case  is  370  ns,  which  is  acceptable  for  driving  a  400-450  ns 
implosion. 


For  scaling  purposes,  resort  to  Massionier's  and  DeMarcu's  analysis 
(Refs.  3  and  4),  suggest  a  cross-sectional  area  for  a  fuse  based  on  the  param¬ 
eters  of  the  driving  circuit  and  on  the  physical  properties  of  the  fuse 
material  of  interest. 


s 


W3/2 

"  V  LJs  M 


(30) 


where 

s  =  cross  section  of  fuse  (m2) 

W  3  stored  energy  (J) 

L  *  total  system  inductance  (H) 

V  =*  charge  voltage  of  capacitor  bank  (V) 

kxa  =  set  of  parameters  describing  the  material;  kxa  =  1.2  x  1017 
(copper) 

For  the  data  In  Figures  16  and  17, 

U  3  200  kJ 
V  3  50  kV 
L  3  67  nH 


Thus,  Equation  30  would  predict 
s  3  7.6  x  10’6m2 

3  0.076  cm2  (31) 

The  fuse  in  question  was  21  cm  wide  and  0.001"  thick  or 
s2  3  0.053  cm2 

or  roughly  70  percent  of  that  predicted  by  Massionier  model. 


3.  Massionier,  et  al . ,  "Rapid  Transfer  of  Magnetic  Energy  by  Means  of 
Exploding  Foil,"  Rev.  Sci .  Instr.  37,  No.  10,  1966. 

4.  DeMarcu,  Berkhart,  "Characteristics  of  a  Magnetic  Energy  Storage  Switches 
Using  Exploding  Foils,"  Rev.  Sci.  Instr.  41,  No.  9,  1970. 
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Scaling  upward  for  a  system  where 

W  =  2  MJ 
V  *  120  kV 
L  =  9.2  nH 

The  predicted  area  is  0.45  cm  or  s2  =  0.32  cm.  For  interest,  a  foil  1  mil 
thick  would  then  be  only  1.3  m  wide. 

Figure  18  shows  a  plot  of  material  resistivity,  p,  versus  specific  energy 
dissipated  in  the  fuse.  The  functional  relationship  between  p  and  specific 
energy  is  open  to  question,  but  for  simple  approximations  we  will  take  the 
empirical  data  of  Figure  18.  Recalling  that  the  previous  analysis  (Fig.  11) 
indicated  that  a  maximum  of  670  kJ  must  be  dissipated  In  the  fuse,  and  taking 
approximately  6  kJ/g  as  the  upper  limit  of  useful  specific  energy  from 
Figure  18,  indicates  that  112  g  of  material  could  be  utilized.  At  a  density 
of  8.9  g/cc  and  a  cross  section  of  0.32  cm2,  this  implies  a  fuse  length 
of  40  cm.  If  the  fuse  reaches  a  maximum  resistivity  of  520  yn  cm,  the  fuse 
that  is  0.32  cm2  x  39  cm  has  a  peak  resistance  of  63  mn.  From  Figure  9,  a 
fuse  with  Rf  of  63  mn  would  drive  an  implosion  to  better  than  400  kJ  of 
kinetic  energy  or  20  percent  overall  kinetic  efficiency.  The  initial  resis¬ 
tivity  of  this  fuse  would  be  0.2  mn,  which  is  much  less  than  the  1  mn  value 
found  to  be  a  realistic  maximum  initial  R...  Note  that  the  interpretation 
attached  to  the  data  in  Figure  18  is  pessimistic  because  the  resistivity 
curve  appears  to  be  clearly  steepening  (not  yet  having  reached  the  plateau 
assumed  in  the  model  of  R^).  On  the  other  hand.  Figure  9  shows  that,  while 
increasing  resistivity  (or  increasing  Rf)  will  help  somewhat,  the  marginal 
gains  are  small. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 

In  conclusion,  simple  extrapolation  of  already  existing  data  leads  to  a 
conceptual  design  for  a  fused  opening  switch  which  can  be  implemented  on  a 
2-MJ  system.  The  resulting  plasma  implosion  should  be  compared  against  that 
which  can  be  obtained  by  directly  driving  the  plasma  from  the  capacitive 
energy  storage.  Using  an  initial  geometry  of  7-cm  radius  and  2-cm  height, 
and  requiring  for  stability  reasons  that  the  direct  driven  implosion  be 
complete  in  less  than  1.4  ys,  results  in  the  coupling  of  approximately  400  kJ 
of  kinetic  energy  to  the  implosion  be  the  direct  case.  This  compares  very 
favorably  with  the  400  kJ  of  kinetic  energy  implied  in  the  above  analysis. 

The  advantage  of  the  inductive  system  is  clearly  the  time  scale  on  which  the 
energy  is  delivered.  The  inductive  system  promises  400-ns  implosions,  a 
factor  of  more  than  3  faster  than  the  1.4  ys  direct-driven  implosions. 
Apparently  significant  gains  can  be  achieved  by  this  fairly  modest  reduction 
in  implosion  time. 


